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A PHASE-MODULATED TRIPLE-EXPOSURE HOLOGRAPHIC INTERFEROMETRIC 
TECHNIQUE FOR HOLOGRAPHIC NON-DESTRUCTIVE TESTING 

Abstract 

A phase-modulated triple-exposure technique is incorporated into the 
holographic non-destructive test (HNDT) system currently existing at NASA 
MSFC. The technique is able to achieve a goal of simultaneously identi- 
fying the zero-order fringe and determining the direction of motion (or 
displacement). Basically, the technique involves the addition of one more 
exposure, during the loading of the tested object, to the conventional 
double-exposure hologram. A phase shifter is added to either the object 
beam or the reference beam during the second and the third exposure. 

Theoretical analysis with the assistance of computer simulation has 
illustrated the feasibility of implementing the phase-modulatioiTand triple 
exposure in the HNDT systems. .'Iain advantages of the technique are the 
enhancement of accuracy in data interpretation and a better determination 
of the nature of the flaws in the tested object. 

Furthermore, the concept of dual hologram fringe subtraction is also 
discussed. The unrealizability of this idea to the HNDT system has been 
demonstrated by results from computer simulation of the subtraction of 
fringe patterns of Fraunhofer diffraction through a single slit and the 
Airy rings. 


A PHASE-MODULATED TRIPLE-EXPOSURE HOLOGRAPHIC INTERFEROMETRIC 

.,! - 

TECHNIQUE FOR HOLOGRAPHIC NON-DESTRUCTIVE TESTING 

I. Introduction 

In holographic non-destructive test (HNOT) applications, It Is the Inher- 
ent nature of the holographic Interference patterns such that a tedlus process 
Is usually required to perform data analysis and Interpretation. In addition 
to the complexity of the Interference fringe patterns which may be given rise 
by the complexity of the test object and the nature of the defects, two other 
problems are often encountered. First Is the Identification of the no-motion 
fringe; second, the detection of the direction of motion. The fringes due to 
no-motion may be considered as a kind of noise that tends to blur the signal 
which leads to the detection of the defects In the test object, therefore, the 
elimination of these fringes Is highly desirable. Direction of motion of the 
object surface Is Important to the determination of the nature of the flaw 
that caused It. 

In order to achieve the purpose of Isclatinq the fringes related to the 

1 

flaws, an Idea was first proposed to subtract one Interference fringe from 
another. Through the subtraction. It was hoped that the residual fringe frag- 
ments, which constitute the difference between any two given fringe patterns, 
will give a simple and clear representation of the flawed regions. Any no- 
motion fringes will be eliminated automatically through the process since they 
are resulted from a common origin. Although this seemed to be an excellent 
idea, after careful studies. It was found that the Idea cannot be realized 
mainly due to two factors: (1) The difficulty In the subtraction of the fringes 
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(2) The Impossibility of the Interpretation of the results quantitatively 
even the subtraction Is assumed to have been carried out perfectly. The 
difficulty In the subtraction of the fringes will be demonstirated by the 
technique of computer simulation as described later In the report. 

Instead of the original Idea of fringe subtraction, we have found that 
a phase-modulated triple-exposure holographic interferometric technique can 
be used to solve the problems of no-motion fringe Identification and the 
detection of the direction of motion simultaneously. Section II will be a 
theoretical discussion of the fringe subtraction and some computer simulated 
results. The phase-modulated triple-exposure holographic technique will be 
provided in Section III. Finally, the Incorporation of the technique In the 
HNDT system will be shown In the discussion part of Section IV. 

II. Fringe Subtraction 

In this section, we shall discuss the concept of fringe subtraction and 
show why it will not be feasible for HNDT applications. 

The basic idea of fringe subtraction Is to subtract one Interferometric 
fringe pattern from another. As an example, one may consider two images, 
one a referen-e set of interferometric fri'nges, and the other a test set of 
••nterferometric fringes. Both of these sets of fringes are obtained by using 
the same object. One may assume that the first set of fringes is obtained 
when the object is in perfect condition and Is considered as the reference 
set. The second set is obtained after the object has been used or Its com- 
position has been modified or degraded. If the first set of the fringes may 
be subtracted from the second set of fringes under ideal conditions a residual 
set of fringes Is obtained. It Is hoped that the residual set of fringes 
resulted from the subtraction can be used to Indicate the change of composition 
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of the object. Unfortunately, it wa''. concluded after careful studies that 
this idea is not practical with regard to its potential applications to a 
HNOT system. We shall present our evidence as follows. 

First we shall investigate three fundamental fringe patterns and dis- 
cuss the subtraction of fringes based on these patterns. 

(1) Fresnel zone plate 

Fresnel zone plates are patterns of concentric circles, as shown in 
Figure 1(a), in which all black and white annuli are of equal area. It was 
found that when two identical zone-plate patterns that are not in concentric 
alignment are superimposed, fringes would appear as shown in Figure 1(b). 

These fringes are normal to the line of relative displacement. As it was 
derived in Reference 1, the number of fringes permilimeter (mm) is directly 
proportional to the displacement: 

d = R^n (1) 

where d is the relative displacement between the two centers, mm; R is the 
radius of the central circle in the zone plate, mm; and n is the number of 
fringes per mm. - 

The fact that moire fringes occur even in the case of a simple zone plate 
subtraction under the condition of misalignment indicated that critical and 
exact fringe alignment is required in our fringe subtraction subsystem. 

However, zone-plate patterns obviously do not represent typical holo- 
graphic interferometric fringes simply because the zone plate patterns consist 
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Figure 1. (a) A Fresnel zone 
plate pattern, (b) Supperposition 
of two non-concentric zone plate 
patterns. 


Figure 2. (a) The Fraunhofer diffraction 
pattern of a single slit, (b) Its 
experimental verification. 


Figure 4. Normalized airy ring 
diffraction intensity pattern. 


Figure 3. Circular aperture geometry 


Figure 5. Airy rings 
with aperture 
diameter equals to 
(a) 0.5 mm; (b) 1.5 mm 
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of only black (opaque) or white (transparent) rings without possessing any 
gray tones. Hence a second case Is considered below. 

(2) Fraunhofer diffraction patterns of a single sllt^. 


The Fraunhofer diffraction pattern (of Eq. (2) below) of a single silt 
Is shown In Figure 2(a) with experimental verification shown In Figure 2(b). 


The normalized Intensity of the pattern may be written as 

2 

me 


1(e) _ r sir 

TW ■ ' 6 


) 


( 2 ) 


where e = (Kb/2) sine 

b = width of the slit 
k = 27 t/x 

A = wavelength of the Illumination 
9 = angle measured from the center of the slit to 

The subtraction of the fringes can be analytically 

2 . , .2 
/ - sine X , . sirUB + AB) x 

{ a ) ( b -J+Js ^ 


the fringe location, 
represented by 

(3) 


where a and b are adjustable constants representing the Intensity scaling factor 
and AB denotes a misalignment factor. 

In order to see the Influence of the parameters a, b, and AB on the fringe 
subtraction, a computer simulation prograni-of the problem Is being Implemented. 

Although the diffraction pattern of a thin slit bears more resemblence to 
that of the holographic Interferometry than the pattern of the zone plate. It 
is not two-dimensional. For this reason, we shall describe a third case as 
follows. 


(3) The Airy pattern^ 

The normalized Intensity distribution of the Airy pattern can be written 


as 



2J^(k a sine) 
— k'd' Sihe 



( 4 ) 
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where Is the zero order Bessel function of the first kind; and 
k » Zk/\ 

\ ■ wavelength of the Illumination .iK- 

a > radius of the circular aperture 
6 - q/R; q and R are shown In Figure 3. 

Equation (4) Is plotted and shown In Figure 4. 

Photographs of the Airy rings of different apertures are shown In 
Figure 5. 

It would also be Interesting to study the difference, S, of the Inten- 
sity of two fringe patterns under the Influence of c, d, and A0: 


- r ri 1(9 

^ ^ TW ■ ° 1(0) * 


(5) 


where c, d are the intensity scaling parameters and A0 Is a factor reflect- 
ing misalignment. 

Four computer programs have been written for the fringe patterns of the 
Fraunhofer diffraction of single slit (corresponding to Equations (2) and (3)) 
and the Airy fringe pattern (corresponding to Equations (4) andjB)). The 
computer program for Equation (2) (sin^x/x^) is listed as program 1; Equation 
(3) is listed as program 2; Equation (4) as. program 3; and Equation (5) as 
program 4. These computer programs are listed ^n Appendix A. 

The function sin^x/x^ versus x is plotted in Figure 5; due to the even 
symmetry of the function, it is only necessary to plot the function for 
negative x. The vertical coordinate of each asterisk indicates the normalized 
intensity of the fringes of the Fraunhofer diffraction. The function of 
fringe subtraction, 

(sinx/x)2 - (A sin(x+d)/(x+d))2 

versus x for positive x are plotted, with various parameters In Figures 6 
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Figure 7 
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Figure 8 
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Figure 11 
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Figure 12 
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through 17. In Figures 6 through 9, the parameter A = 0.5, and d ■ O.lir, 

0.27T, 0,3ir, and 0.5ir respectively. For the same series of value of the 
parameter d. Figures 10 through 13 are with A = 0.8; and Figures 14 through 
17 are with A = 1.0. Figures 6 through 17 clearly demonstrated that slight 
misalignment and/or misadjustment of the intensity of two simple fringe 
patterns will give rise to another fringe pattern which is just as compli- 
cated irf relative fringe contrast as the original fringe pattern. The 
computer-simulated result implies that when the holographic interference 
fringe patterns of real (’oubl e-exposure holograms or real-time holograms 
are combined without any precise knowledge of the phase and amplitude para- 
meters, uninterpretable and confusing result will easily occur. 

The outputs of computer programs 3 and 4 showing the fringe subtractions 
of the Airy ring fringe patterns demonstrated the same fact; further rein- 
forced the correctness that the idea of fringe subtraction is not practical 
for the HNDT work. 

Based on the above conclusion, a new approach was searched for the 
purpose of achieving the original goal to improve the data reduction process 
in a HNDT system. A technique called the phase-modulated triple-exposure 
holographic interferometric method was found. This technique will be discussed 
in the next section. 

III. Phase-modulated Triple-exposure Holographic Interferometry 

3 

The technique is based on D. Gabor ^^'s suggestion that complex 
addition of images or the interference patterns is a type of "true" optical 
filtering process. Complex image addition or subtraction may be called an 
"optical image synthesis" process. We shall utilize this concept together 
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with a phase shift operation during a triple-exposure process to identify 
the no-motion fringe; in the meantime, it is expected that the direction 
of the object's displacement may also be determined. ..jv- 

The intensity of the fringes in holographic interferometry is propor- 
tional to the square of the characteristic function. In general, It Is nec 
essary to solve the equation 


Kr) = 





exp [j ^ (r, t)] dt 


( 6 ) 


\ '0 ^ 

where IqCt) denotes the intensity at point r as reconstructed from an ordinary 
single-exposure hologram, T, the exposure time, and ^ (r, t) the phase of the 
object beam relative to the reference beam. 

We shall assume that in general, the loading of the object causes the 
phase (j> (r, t) to vary according to 


♦(r, t) = 


♦i(r) ; t < t^ , 

f(r, t) ; < t < tg , 


(7) 


Or. 

U. 

♦ (r, t) = (j.^(r) [u(t) - u(t - tg)] 

+ f(r, t) [u(t - tg) - u(t - tg)] 

+ ♦2^'"^ 

where u(t) is the unit-step function defined by , 

u(t) = 1 , t ^ 0 

= 0 , t < 0 . (9) 

The phase function (^(r, t) is plotted and shown in Figure 18. In Figure 
18, the broadened parts of the curve show the time intervals while exposures 
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of the hologram are made. These time intervals are between t^ and t^; t^ and 
t^; and t^ and tg. Ordinary double*exposure holograms are made usually between 
t^ and t^ ; and t^ and t^. A third exposure between t^ and'k" is added in the 
present case. This additional exposure during the loading process is the key 
to the triple exposure technique. 

According to the exposure scheme shown in Fig. 18, Eq. (6) may be written 
as 


I(r) 


lot--) 


f 


V 



exp[j(|>j(r)] dt 



exp[jt(r,t)] dt 



where if we define 



^8 ■ ^7 “ ^3 * 

then in Eq. (10) 


( 10 ) 


h 


( 11 ) 


T = Tj + * T 3 (12) 

In addition to the third exposure, we may also introduce various constant 
phase factors by inserting phase-plates into the laser beam during the second 
and their exposures. If these phase factors are denoted by a and 8 respectively. 
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then Eq. (10) becomes 



The aoove equation represents the basic principle for the phase-modulated 
triple-exposure holographic interferometry. 

To illustrate the effect of Eq. (13), we shall present an example below. 

The function f(r,t) in Eqs. (7) and (8) is generally unknown^ and position 
dependent, and therefore is subject to one's conjucture. We shall assume, for 
simplicity, that it is space-independent and varies linearly with respect to 
time. In addition, we shall assume that the exposure during the loading process 
begins while the loading begins/and ends at the end of the loacffng period. Con- 
sequently, the phase function becomes 

♦(r); t'^’t 2 , 

♦ (r,t) = < [(t - tj)/(t 2 - tj)] A<|.(r) ; t^ < t < tg , 

(<i(r) + A(|)(r) ; t > tg . (14) 

Equation (14) is plotted and shown in Fig. 19. The three exposures are 
now between time periods t^, t^; t^, tg; and tg, t^. 

For simplicity, we shall let 


tj - t, - T, , 

^3 * ^2 * '2 * 


‘c - ‘3 • T. 
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Phase function (j>(r,t) versus t for the 



where T ■ 2Tj + Tj 

By letting a > B ■ 0 and substituting Eq. (14) Into Eq. (13)> one obtalns^'^ 
I(r) - [Ij,(r)/T2] [2Ti cos X + Tj (sin X/X)]2 (16) 

If Instead a « B ■ 0, we have a > 0 but B ■ w, Eq. (16) should be 

I(r) - [Io(r)/T2] [2Tj sin X)* + (T^ sin X/X)2] ; (17) 

and for a ir/2, B ■ w, and -A4, Eq. (16) becomes 

I(r) = [I^,(r)/T2] [2Tj sin x - (sin X/X)]2. (18) 

The variable X In Eqs. (16) through (18) Is defined by 

X = ^ . (19) 

Equations (17) and (18) can be normalized respectively that yields 


respectively 





T(r) = sin2 x + sln^ X/X^ , 


(20) 

and 





T(r) = [sin X - D sin X/X]2 , 


(21) 

where 

• 




T(r) = T2l(r)/[4Tj2 I^(r)] , 


(22) 

and 


r- 



D = V(2Tj) 


(23) 


U 


Computer programs have been written to calculate and plot Eqs. (20) and 
(21). Corresponding to Eq. (20), Figures 20 through 23 show the variation 
of the normalized Intensity In the reconstructed Image against the phase 
change X = a^/2 for 0 = 1, 2, 3, and 4 respectively. It Is clear that all 
dark fringes occur at an Interval of 2it for all values of 0. The zero-order 
fringe is reconstructed much brighter than higher order fringes as shown in 
Figures 21 through 23 for 0^2. Hence, it may be concluded that the zero- 
order fringe due to no-motion can be Identified by having a change of phase 
of IT In one of the beam for the final exposure and a long Intermediate 
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Figure 21 
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Figure 22 
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exposure during the loading process such that D » T2/(2Tj) ^ 2 . This 
technique, however, will not detect the direction of motion. The detec- 
tion of the direction of motion is discussed below. 

Corresponding to Eq. (21), Figs. 24 through 27 show the normalized 
intensity for D * 1, 2, 3, and 4 respectively versus the phase change X. 
Again, it may be noted that all dark fringes occur at an interval of 2ir 
for all D, and the zero-order fringe is much more intense compared with 
the rest. For the case where Ai|) is negative, the intensity of bright 
fringes gradually increases with the increase of the fringe order. An 
increase of 3 times from the first-order fringe to the 5th-order fringe 
is visible for D = 1; and about 2 times for D = 2. For the case where ^ 
is positive (plots not shown), the intensity of bright fringes decreases 
with the increase of the fringe order. Since the sign of the phase A4 
is correlated to the direction of motion, this characteristic of the 
fringes enables one to discern the positive direction of motion from the 
negative direction of motion of the object at position r. 

w. 

t 

The computer programs given the above plots are listed as program 5 
and program 6 in Appendix A. ' — 
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IV. Discussion 

From the examples given in the last section, it can be seen that the 
examples in Sec. Ill indicate that the phase-modulated triplifi-exposure holo- 
graphic interferometric technique can be used in the identification of the 
zero-order fringe and the determination of the direction of the displace- 
ment. Although the phase variation during the loading process was assumed 
to be linear with respect to time, the result would not be qualitatively 
different if other time dependence was assumed. 

The incorporation of this technique into a HNDT system can be made by 
applying a removable phase-shifter to the object beam. A block diagram 
showing the phase-modulated HNDT system is given in Fig. 28. The operational 
procedure is described below. During the first exposure, the phase shifter P 
is removed from the object beam. The phase shifter is replaced into the 
object beam during the second and third exposures. The amount of phase shift 
may be different. For example, the phase shift during the second exposure 
may be it and that during the third exposure may be tt/ 2. It should be noted 
that the phase-shifter may also be inserted into the reference beam instead. 

If a flat (or quasi-flat) object isjjnder test, the analysis of data can 
be made by concentrating our attention on the higher order fringes and use 
the principle that the fringe pattern is of the following form^ 

A{x^ + y2) + Bx + Cy = (n - -j)x , (24) 

where \ is the wavelength of light; n is a positive integer greater than or 
equal to 1; A, B, and C are system dependent parameters. It may be concluded 
that due to the fact that the zero-order fringe is not counted, a more 
accurate assessment of surface displacement in the HNDT system can be achieved. 
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Figure 28. A block diagram for 
the phase-mo.-'ulated HNDT system. 
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APPENDIX A 


Computer Programs 



1 

2 

3 

5 

6 

7 

8 
9 

10 

11 

12 

13 

U 

15 

16 
17 
IP 

19 

20 
21 
22 

23 

24 

25 

26 
27 


CPROGRAM 1 

DIMENSION P(1C2),F<100) 

INTEGER PtBLAKKtSTARtDOTtDYtZERO 
DATA PLANK, STAR, DOT/ ^ 

DATA ZERO/'O'/ 

FM1N«99*>,9 
F«AX«-999.9 
DO 10 1=1,41 
DX»(X-1)«7.14159«0.1 
F (1)«(SINCDX) /DX)**2 
IF(I.EQ.I) F(l)=1.0 
1F(F(I) .GT. FPAX) FMAX«F(1) 

IF (F(l) .LT. FMIN) FMIN=F(I) 

10 CONTINUE 

P <1 )* DOT 
P (102) = D0T 
DO 20 1=2,101 
20 P(I)=eLANK 

V.RITE(6,ir)0> FriN,FMAX , 

100 FORPATC '1',1CX,'PL0T OF SINX**2/X*«2 VERSUS X't/ ' ',5X ,‘'S1NX = *2 
C/X«*2',F6.2,95x,FE.2) 

SCALC»(FMAX-FMN)*C.001 

MT=(-FMIN'/SCALE*1 

DY=50 

F1»FMIN'»(50^MY)*SCALE 
FJ«FMIN*(i»Y-5C)»SCALE 
P (DY )*ZERO 


28* 


V RITE (6, 104) F J, (P (J ),J*1, ICC) f FI 

29* 

104 

FOBMATC' ',16X,F5 .3,100A1,F5.3> 

30* 


P(DY)*0LANK 

31* 


WRITE (6,101) 

32* 

101 

FORKAT {'0',1CX,115('-'); 

33* 


DO 3u 1-1,81 

34* 


DX=(I-1)*3. 14159*0. 05 

35* 


DY=(F(I)-FMIN)/SCALE*1 

36* 


DY=DY-MY*50 

37* 


IF(DY.GT.IOI) GO T0'“j40 

38* 


if(dy.lt.i) gc to 40 

39* 


P (0 V )»STAP 

40 • 

40 

WR1TE(6,102) DX,(P(J),J-1,102) 

41* 

102 

FORMAT ( ' ',5X,F12.5,4X,1C2A1) 

42* 


IF(DY.GT. 101)60 TO 30 

43* 


IF (0Y.EO.1) r(DY)*DOT 

44* 


IF(DY.GT.I) P(DY)«BLANK 

45* 

30 

CONTINUE 

46* 


WRITE (6,101) 

47* 


WRITE(6t103) 

48* 

103 

FORMAT ('1') 

49* 


STOP 

50* 


FNN 




_ 


C program 2 


2* 


DIMENSION P(1C2).F(80).A(3),0(4) 

3« 


INTEGER P. blank, star, DOT, 0V«2eR0 

A* 


DATA BLANK, STAR, dot/ ' 

5* 


DATA 2ERO/'0'/ 

6* 


DATA (A(I) ,1*1,3)/ 0,5, 0.8, 1,0/ 

7« 


DATA (D(l) .1*1,4) /C. 01,0, 02, c.:.3,o. 

8* 


EMIN-999.9 

9* 


F’1AX*-999.9 

10* 


DO 5 1A*1,3 

11* 


DO 7 1J»1,4 

12* 


DO 10 1*1 ,61 

13* 


DX*(1-1)*X.U159*0.1 

U* 


DJ*D (IJ)*3. 14159 

15* 


B*(A(IA)*S1N(DX*DJ )/(DX*DJ))**2 

16* 


MI)*(S1N(DX)/DX)**2>B 

17* 


IF(I.EQ.I) F(I)=1-B 

18* 


'IF(Fd) ,6T. FMAX) FHAX*F(1) 

19* 


IF (F(I) .-.r, FHIN) FMIN*F(X) 

20* 

10 

CONTINUE 

21* 


P(1)» DOT 

22* 


P (102)=DOT 

23* 


DO 20 1*2,101 

24* 

20 

P(1)*BLANIC 

25* 


WRITE (6 ,100) 

26* 

100 

FORMAT ('1 ',1 CX.'PLOT OF (SINX/X)** 

27* 


ex') 

2£* 


WRITE (6, 99) A(7A),0(IJ),FHIN,FMAX 

29* 

99 

FORMAT (' 1CX,'A=',F5.2,10X,'D*' 

30* 


C ,F10.4,3X,'F7AX=',F10.4) 

31* 


SCALE=(FMAX-FPIN)*C,0C1 

32* 


MY*(-FMIN)/SCALE*1 

33* 


DY*50 

3A* , 


F J»FMIN*(MY-5C)*SCALE 

35* 


FI=FM1N*(50*MY)*SCALE 

36* 


P(DY)=ZERO 

37* 


WRITE (6, 104) FJ,(P(J),J*1,10C) ,F1 

31: * 

104 

FORMAT!' ',16X,F5.3,100a1 ,FS,3) 

39* 


P (dy)=blank 

40* 


WRITE (6,101) 

41* 

101 

FORMAT ('C',1CX,115('-'')) 

42* 


DO 30 1 = 1 ,41 

43* 


OX=(I-1)*3. 14159*0.1 

44* 


DY=(F(I)-FHIN)/SCALE*1 

45* 


DY=D Y-HY*50 

46* 


IF(DY.GT.IOI) GO TO 40-j 

47* 


IF(OY.LT.I) 6C TO 40 ' 

48* 


P(DY)=STAR 

49* 


P(DY)=STAR 

50* 

40 

WR1TE(6,102) DX,(P(J),J*1,102) 

51* 

102 

format ( ' ',5X,F12.5,4X,1C2A1) 

52* 


IFCOY.GT. 101)60 TO 30 

53* 


IF (DY.E0.1) P(0Y)*D0T 

54* 


IF(DY.GT.I) P(DY)*BLANK 

j5* 

30 

CONTINUE 

56* 


WRITE (6,101) 

57* 


WRITE (6,103) 

58* 

103 

FORMAT ('1') 

59* 


FMAX*-999.9 

60* 


FM1N»999.9 

61* 

7 

CONTINUE 

62* 

5 

CONTINUE 

63* 


STOP 

64* 


END 


"1 
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I 


Ui 

n < *4 


FROGHAM No. 3 


.i| . 


1 * 
•5 * 


6 * 
7 « 


c 

r 

•II', 


C PROG PA'* 7 

C lr<?.\s inrr PCI “Z ) , f 
I\'TEGER F.SLA '.K , £ ta r, pot, d y 

_ data ELA^'',STAF,&r7/ ' 

F'''Ifi = 9/0 ,0 ^ 

F--*AX = -9 39. 9 

DO ! = ■; ,51 

CY = Cl-1).nr4 

F0=°3SL(:^,‘!) 


'l 1 * 
1 
1 

in 

IF(FCI) ,0T, F.VAX) FKAX=i 
IF Cr(l) ,LT. F •'!*;) F?‘I\; 
C n.NT I \'U = 

FCI) 

= FCI) 

1 ^ * 


P (1 ) = DOT 

- — 

1 ' * 
-iF • 

■i 7 » 

- 

P Cl 02 ) = 00T 
5 '' 2 j I = 2 , 1 •: 1 


i ' ^ 

^ 1 
4 . . 

P Cl )=r-LA.\‘' ~ 


.1 = * 

10G 

U'^ITE (6 , irQ) fmIn.FMAX 
F'^KMATC 'I'.I'x.'plOT OF 

S I \ X* 

^ IT 


C/x**2',Fc.Z,95x,F2.P'> ~ 

•£1 * 

2?* 

101 

i'p:t£ C6,'^.j1) 

_f OP.VAT C'-', 1 •'X, 1 '! 5 J 


Z 7 ^ 


s C A L 2 = ( F ;< A X - F : N )■* C . C 1 ' 


2 ^ * 
*> c ^ 

L. * 

> < * 


D'' 3 j 1 = 1 ,51 
DX=CI-1 ) *'•.4 




DY=CFCl)-Ff'IN)/scsLE-*-1 

“ 

^ ^ IT 

Z'^ * 


P cdy)=stao 

C6 ,1 '•2) PX,CPCJ),J=1 

, 1C2) 

20- 

102 

FORMAT C ' :X , F1Z.5 ,4X, 

lOZll 

3 ** * 


I ■= CdY.FG.1) fc.dy)=pot 

■ 

j » 


IFC0Y.GT.1) ? CD Y) =CLANK 


i?*' 

30 

C 0 \ T i n u e 



VERSUS X',/ 


»5x,'sinx**2 


:>o- 


W ^ I T = ( 6 , c 1 ) 
URIT£(a, 1 'i 3 ) 
FORMAT C'^') 
STOP 
EN’O 


D- CF COMPILATION : 


NO diagnostics. 


ORIGINAL 

D/wm OLaL-Ixii 



PHOGHAX No. 4 


’ * 
•' . 

* « 

4 • 

5 . 

C 

P«0 OR M 

OI.VFNSION PC 
I T r u £ R F , L 
u\T f 1 L A N , S 

O.ATA (AC), I 

• 5 • 


oata C0CI),T 

'• • 


Ff INrOCP, C 

c • 


, FVAX=-999.9 

50 ' IA=1 ,3 

' r , 


1''. 


50 Y u = i ,4 

1 1 • 


DO 1 J 1=1.51 

K- 


7X=(I-1 .")* 

1 !• 

U« 


FS=PSSLCDY,31 
0J=0CJ) 

tx=cx*o J 

1 '• 


1 -f * 


"^'0 = rSSL CTX', 

1 7 » 


r = (?-A(iA)'r.j 

I'*- 


r Cl ) = (2*c*‘/?> 

« A 
1 ' ' 


If(cCl) .GT. 



IF CfC) ,LT. 

CO.\Tlfjjp 

21* 

1C 

:?• 


PC1)= DOT 

22' 


P Cl '■'2 ) = c C’’ 

^ 4 » 


00 2 2 I = 2,1 PI 

2 ^ « 

2C 

P Cl )=£.LA\k' 


- • 

jC* 
5 * ' 

jC •• 
i 3 • 
3 i « 
5*^- 
i ^ • 
3'’ ■ 




..5 * 
4.5 • 

.,7. 

hT* 

5-'« 
51 • 

5 r * 


.) 

'/ 


;-r 




CO 


101 


132 

30 

103 


IT£ (6. IOC.) _ 

^F^^KrAT 1 -X.'PLOT r.F (SI.\X/ 

w ^ I T £ f 6 , 0 ) A ( I A ) , D c IJ ) , F>: I \ . P .V, 
^0^.':AT (- ^^CX,-A=^F5.^.1JX. 
C •f"‘’j.H,3'','F'AX = ',F10.A) 

J'-ITE (£,‘'C1) 

Format ('"',1 'X.IISf'-')) 
SCALt=.(F'iAX-F*-If.-)*C.C1 

*?0 1=1,5 1 

uX=(:-1).^,4 

DX=(F(I)-c,vi.v)/sf4LH^1 * 

r Coy ) = sTA7 

w^iT£ (‘ , 1^2) ry, vp(j ) .jrT , i'52) 
format C ' '| 'X , F 1 4. 5 ,4X, 1J2 A1 ) 

-At COY.cw.1) fCOYlsOOT^ 

I FCf'Y.GT. '• ) PCOY)=tLANX 
C ''MI VUE 

w?IT£ (6,101)- 
W‘^IT£(^,1':5) 

format 

FVAX=-999.c_ 

F'* I\ = C90 ,c 
C ONTIVUE 

CONTINUE , 

STOP 
END 


x)** 2 -(, 


A X 
'J = 


t F5 


OF COMP I L A, i ON: 


NO DIA gnostics. 


4l 





I--- 


:>* 

4 * 

» 

b* 

/* 

a* 

iu* 

U*- 

i:>* 

Lh* 

lb* 

io* 

17*^ 

io» 

«:u* 

<ii» 

2 b* 

<: 4 » 


C l^r%Oui\Ai*i b 

Lil.*iLNblON H(iU2) »FUbO> »D(4’ 

IimIcocK P»oi_AuK»blMR»DOT»L)Y 
DmiA dLANK# STaK »U0 I / • •»*♦»»».*/ 

UAfA ( 0 ( 1 ) »I = l» 4 )/i, 0 » 2 . 0 »b. 0 » 4 . 0 / 

PiXl.\=99y.9 
FMmX=-999.9 
Do 7 U=1»H 
DO 10 1=1»121 
UX=( I-l)*b.l4io9*u.l 
a=lJ**2*(SlN(DX) )**2/(DX**«l) 

K(i) = (SIfJiOX) )**24 u 
I^(I*L 0«'l) F(i)— 2«u 
ir^Fd) .ol. F.V.AX) FMAX=F(i) 

IF (F(I) .lT. FMifM) FMIN=F(I) 

lo COi^iTInUE 

SCmLE=(FMaX-FM1N)*u.01 
P(i)= DOT 
Plx02)=DOr 
Du 20 I=2»101 

20 P(j.)=tJLANK 

FrtMX=100*bCALt 
tiiU.iTE(6» luu) 

iou FUKX.AT ( *1* rl(JX» VuOr OF ( SlivIX) ♦♦2+D**2 (SINX) **2/X**2 VERbUS X •) 
«»Ki.Tl.( 6»99) D( IJ) rr*'.; nj»FMAa 

99 FoKNiAT (• * »1UA» »F5.2» • *»/» * »14X» »FMIN=» »F10.4»bX» »FKAX» 

C»FiU.4) 

WRiTE(6*lUl) 

101 FOKiMAT ( * 0 • » Ibx » ilo ( * “• ) ) 

Uu oO I=i»l21 
OX-(I-l)*b.l4ibQ*0.1 
DY= ( F ( I ) -FMIN) /SCAL.E41 
IF (DY • GT« lUl ) bO Tb 40 
P(LiY)=STAK 

40 WKiTE ( 5» 102) Da » ( p ( J ) » J=1 » i02)~ 

102 FO.XF.AT ( * • »bX»Fl^.5»4X»102Al) 

IF (DY.EQ.l) P(DY)=00T 
IF(DY.GT.101)GO Tu 3U 
IF(DY.GT.l) P(uY)=pLANK 

30 Continue 

mKITE (o»1o 1) 
kVRiTE (6r lU3) 

105 FuRMAT (»j.*) nrnnixTAT t>app m 

FHAX~* 999»9 ORIGDTAIi PAGE lo 

FM 1 N= 999 . 9 OF POOR QUALITY 

7 CuRTiNUE 

STuP 


CGwPiLAv ion: 


DIA6i>»0bTICS. 


1 



1 + 

C PuOoUAi't D 



Liii'iENb 10\ H ( lo2) » K <2o0 ) » J ( ‘t ) 



Irt 1 cbc-H P » oi_Ai‘iK » o 1 mR » DOT » OY 

»«» 


OaIa u^-AN^»bTMR»oui/ * •»*♦»»»• 



UAFA (0(1) »I=i»4)/x*0»2»0»«i«0»4 

u* 


Ft*‘lu=99b • 9 

/* 


Fi'*’.AX=-999« V 

ij* 


OU 7 U=l/4 

SJ* 


JO lU l=i»241 

io* 


0X-(I-13)*O.14ib9*o.l 

!!♦ 


o=iJ*SlSljX)/uA 



IF (I»cQ«l) O=i0 

io» 


r (i) = (SIM(l.X)-d)**^ 

i*t» 


lF(F(i) .OT. FViAa) F4AX=F(I) 

ia'*- 


IF (F(l) .LT. F.vii>4) FV.IN=Fd) 

Xo* 

iu 

CO.mTINUE 

1/ + 


^OrtLL— ( FMaa“FisIN) ) ♦u • 0 1 

io* 


PiX)= COT 



P(iU2)=D0T 

eLO* 


jo 2U I=2ri01 


Ho 

P(x)=bLANK 



Fi^uiX— lOO*bCALL 



ki.*<i V k£ ( 6 r loo ) 


loo 

FoKiXmT (*i*»luA*' FLOT of (SIimX' 

C.O* 


;.Rx rE(o»99) D( IJ) »FMxM»FMAa 

^o* 

Vb 

F0i\i*iA7 (* ' » lUA» ' j— * » Fb«2» * *», 

<L'l* 


C»FiJ.4) 

<ia* 


.'.t\ iTtl ( vS » IGi ) 


lul 

Fv^K.-.hT ( ' 0 ' » luA » iio ( * “ » ) ) 

Ou ♦ 


DO iU I=lr 24 l 

Ox*^ 


OA’=(I-l)*0.141ijP + o»l 

Oc.* 


JY=(F(i)-F,>.iTN)/So„ut + l 



IF i l»Y • GT • X 0 1 ) oo i j 4ij 



P(oY)=STAR ' 

6o* 

ho 

.-..<1 r j(o» loc) jx» (plj) » J=lri0'2) 

Oo^ 

x02 

ruiOlAT ( * * tbX»Fi4i.b»4X»i02Al) 

J/* 


IF (jY.EQ.l) P(DY)-D0T 



iF(uY»GT«Xol)si ;0 To 30 

Ob* 


IFljY.GT.x) P(jY)=oLANK 

4-U* 

00 

CwuT I^J£ 

4i^ 


iiVKxTE (o»ioi) 



Xr<iTc(o» lUo) 

‘to* 

ioo 

F ok.M A T ( • i • ) 

•ts^* 


FMmX=-999« 9 

4o* 


Fi*ilN=999» 9 


V 

Co»i\jT i^^UE 

4/* 


bToP 

4o* 


Li\u 


U or COf-.PlLAi iOi\; ,\0 l^ImO.<iOST1CS. 


SInX/X )*♦?. VEHbUf; X •) 

• » 14X» 'FMIN=» »H0.4»3X» 'FMAX* 
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